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Abstract. The Doppler effect, the shift in the frequency of sound due to motion, 
is present in both classical gases and quantum superfluids. Here, we perform an in- 
sit% minimally destructive measurement, of the persistent current in a ring-shaped, 
superfluid Bose-Einstein condensate using the Doppler effect. Phonon modes generated 
in this condensate have their frequencies Doppler shifted by a persistent current. This 
frequency shift will cause a standing-wave phonon mode to be “dragged” along with 
the persistent current. By measuring this precession, one can extract the background 
flow velocity. This technique will find utility in experiments where the winding number 
is important, such as in emerging ‘atomtronic’ devices. 
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1. Introduction 

Ring-shaped Bose-Einstein condensates (BECs) use topology to exploit one of the 
key features of a BEC: superfluidity. In particular, the topology supports superfluid 
persistent currents [1], As a result, a ring-shaped condensate forms the basis of 
several so-called ‘atomtronic’ devices: simple circuits that resemble counterparts in 
electronics [2, 3, 4, 5, 6]. The addition of one or more rotating perturbations or weak 
links into the ring can form devices that are similar to the rf-superconducting quantum 
interference device (SQUID) [2, 3, 4] and dc SQUID [5, 6]. Operation of these devices 
typically requires measuring the persistent current. Here, we present a technique for 
measuring the persistent current of a ring that uses the Doppler effect and, unlike other 
methods, is done in-situ and can be minimally destructive. 

Superfluids can be described using a macroscopic wavefunction (or order parameter) 
xp = where n is the density of atoms and f is phase of the wavefunction. In this 

picture, the flow velocity is given by v = (Ti/m)V4>, where m is the mass of the atoms 
and % is the reduced Planck’s constant. Because the wavefunction must be single valued, 
the integral f • dl must equal 27r^, where £ is an integer called the winding number. 
This quantization of the winding number forces flows around a ring of radius R to be 
quantized with an angular velocity Dq = fi/mBf. Thus, the angular flow velocity must 
satisfy Vt — £Dq. 

In addition to supporting bulk persistent flows, the BEC also serves as a medium 
in which sound can travel. Because of the trap’s boundary conditions, only certain 
wavelengths of sound are permitted. These phonon, or Bogoluibov, modes are the 
lowest energy collective excitations of the condensate [7, 8]. In this work, we excite a 
standing wave mode with wavelength A = 27rR, as shown in Fig. l(a-b). Such a standing 
wave is an equal superposition of clockwise and counterclockwise traveling waves with 
the same wavelength. In the presence of a background flow, the Doppler effect shifts 
the relative frequencies of the two traveling waves. This shift causes the standing wave 
to process, as shown in Fig. 1(c). Here, we use this precession to detect the background 
flow velocity of the superfluid. This method is analogous to earlier techniques, where the 
precession of a quadrupole oscillation was used to measure [9, 10] the sign and charge of 
a quantized vortex in a simply connected, harmonic trap, as suggested in Refs. [11, 12]. 

In a ring-shaped condensate, two other methods have been used to measure the 
background flow velocity by determine the winding number £. Both are inherently 
destructive because they require that the BEC be released from the trap. First, sufficient 
expansion of the condensate can yield a hole at the center of the cloud whose size is 
quantized according to £ [13, 14, 15, 16]. Second, experiments releasing both a ring and 
a reference condensate produce spiral interference patterns that indicate £ [17, 18]. 

Our method of detecting rotation by the Doppler effect is sufficiently precise to 
distinguish adjacent values of £. However, it is minimally disruptive: it requires only 
exciting a sound wave and imaging the resulting density modulation. By imaging the 
density modulation using a minimally destructive imaging method (such as dark-field 
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Figure 1. (a) Image of the ring condensate in-situ without an applied 

perturbation, (b) A sinusoidal perturbation excites a standing-wave superposition of 
counterpropogating phonon modes. The white line intersects the ring at the maximum 
and minimum of the perturbation. The left image shows the density; the average 
density is subtracted from the right image. These images show the resulting density 
modulation, 100 /is after the perturbation has been removed, (c) Full density (top) 
and average-subtracted-density (bottom) images taken 9.1 ms after removal of the 
perturbation. The density modulation rotates relative to the initial perturbation in 
the presence of superfluid flow. The winding numbers are shown below the images. 


dispersion imaging [19], phase contrast imaging [20], diffraction-contrast imaging [21], 
partial transfer absorption imaging [22], or Faraday imaging [23]) one can create a 
winding number measurement that is also minimally destructive. Such measurements 
could allow easier implementation of experiments where rings are stirred by weak links 
more than once, such as hysteresis experiments. By measuring the winding number 
in a minimally destructive way, the quantum state of the ring can be known after 
each stage of stirring. One can use multiple minimally-destrnctive winding number 
measurements to increase the sensitivity of a rotation sensor based on the atomtronic 
rf-SQUID. Our system [4] has a sensitivity of 0.1 Hz for a measurement done once 
every 30 s, which includes creating the condensate, stirring, and measuring the 
winding number. This leads to an effective sensitivity to rotation of r\j 5 Hz/yni. 
By repeating stirring and winding number measurement on a single condensate using 
a minimally-destrnctive technique (which yields the same information as a destructive 
measurement), the sensitivity improves by \/Vn^, where V^eas ~ 10 is an achievable 
number of minimally-destructive measurements. 

Finally, it is possible to use this technique to detect rotations of the inertial frame 
of the condensate itself, using the Sagnac effect):. Ref. [24] measured the noise level 
associated with such an excitation-based rotation sensor, finding a sensitivity in their 
system of roughly 1 (rad/s)/yHz. However, their experiment was not configured to 
produce rotation, so the effect was not measured experimentally. Here, we demonstrate 
the feasibility of this idea, by detecting the quantized rotation of the superfluid in the 

f We note that unlike the ring-shaped condensate considered here, not every mode in a simply 
connected condensate is sensitive to both a rotating frame and to vorticity. In particular, the dipole 
oscillation mode is sensitive to a rotating frame but not sensitive to vorticity [11]. 
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2. Theory 


To understand the Doppler effect and its effect on phonons in the ring, let us first 
consider a standing-wave phonon mode in the absence of a persistent current. To 
simplify the calculation, let us consider the equivalent problem of a ID gas of length 
2TrR, where R is the radius of the ring, and let us impose periodic boundary conditions. 
A standing wave mode can be generated by a perturbation of the form V = Vq cos{qO), 
where g is a positive integer, Vq is the amplitude of the perturbation, and 6 is the 
azimuthal angle. This perturbation generates a modulation in the density of the 
form ni^{9) = Ueq -|- 5ncos{q0), where Ueq is the equilibrium density without the 
perturbation and 5n is the amplitude of the density modulation. Upon sudden removal 
of the perturbation, the density is projected onto the spectrum of phonon 

modes. The initial density modulation is then described by the superposition of two 
counterpropogating modes cos{u){q±)t — q±0). Here, indicates a mode traveling 
counterclockwise and g_ indicates a mode traveling clockwise. In the absence of a 
persistent current the two modes have equal frequency: Oi>{q±) = oJq. Consequently, the 
density modulation will oscillate in time, Sn{t) oc cos(a;ot) cos(g6*), without exhibiting 
precession. 

The presence of the current will remove the degeneracy between the two modes 
through the Doppler effect. The frequency of the two modes, in the presence of a 
superflow of velocity v = ifloR, will be given by 

oj{q±)=Uo±^v. ( 1 ) 

The density modulation is then proportional to 


cos[g+0 — uj{q^)t] + co8[q-9 + a;(g_)t] 


= cos 




cos(a;ot) 


( 2 ) 


and hence evolves in space. The azimuthal location of an anti-node of this density 
standing wave precesses around the ring as 
dfrn 


dt 


mR? 


( 3 ) 


Measurement of 4>ra{l) would consequently provide a direct measurement of L We note 
that dfrn/dt does not depend on wq; therefore, the speed of sound and the details of the 
phonon dispersion curve are not relevant in predicting the precession. 

While the above arguments may seem to apply only to a ring, they can easily be 
generalized to any quasi-one-dimensional geometry. In such a case, the shifts in the 
phonon frequencies will still be given by Eq. 1, but with a v that may or may not 
be quantized. For an infinite one-dimensional channel, for example, the flow is not 
quantized. However, two equal but oppositely directed phonon wavepackets (similar to 
those generated in Ref. [25]) traveling in the channel will move with different velocities 
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Figure 2. (a) Schematic of the setup used to trap the atoms and create the double 

ring potential and sinudoidal perturbation. A DMD is illuminated with blue-detuned 
light and is imaged onto the atoms using a telescope (solid lines) , which ensures that 
incoming collimated beams remain collimated after they emerge from the telescope 
(dashed lines), (b) The halftoned patterns || written to the DMD for the bare ring 
potential (left) and the perturbed ring potential (right), (c) After convolution with 
the point-spread function of the imaging system, the potentials formed are smooth, 
and show a clear sinusoidal perturbation (right). 


in the presence of a background flow. The use of a ring allows for a straightforward 
means of detecting the frequency shift: precession of a phonon standing-wave mode. 

3. Experimental details 

We create a ^^Na condensate in a crossed-optical dipole trap, as shown in Fig. 2a. A 
blue-detuned double sheet beam, formed by focusing a TEMqi mode tightly along the 
vertical direction, provides vertical confinement. Confinement in the horizontal plane is 
generated using another blue detuned beam that is shaped using a Texas Instruments 
LightCrafter 3000 digital micromirror device (DMD) §. We position the DMD in our 
imaging system to directly image the surface of the DMD onto the atoms. (This method 
is similar to the photomask method used in Ref. [26], with the DMD replacing the 
photomask.) 

The DMD generates a double ring trap. These experiments use only the inner ring, 
which is shown in Fig. 1(a). In general, our condensates contain 7 x 10® atoms, with 
approximately 5 x 10® in the outer ring and 2 x 10® in the inner ring. The outer ring, with 
mean radius 31(1) /rm || , can be used as a phase reference to measure the winding number 
of the inner ring [27]. The inner ring has a chemical potential of pi/h = 3.1(6) kHz. The 
vertical trapping frequency is 1020(30) Hz and the radial trapping frequency in the inner 
ring is 310(10) Hz. 

In order to generate a persistent current in our ring, we apply a rotating potential 

§ The identification of commercial products is for information only and does not imply recommendation 
or endorsement by the National Institute of Standards and Technology. 

II Unless stated otherwise, uncertainties represent the la combination of statistical and systematic 


errors. 
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generated by a blue-detuned beam steered by an acoustic-optic deflector (AOD). A 
previous paper [4] describes this stirring technique. We verify that this produces 
the desired i state 95% of the time by performing a fully-destructive interference 
measurement [17]. The phase winding in the inner ring can be measured by interfering 
the inner ring with the outer ring condensate, which serves as a phase reference. Fig. 2(d) 
shows a typical interference pattern, which indicates i = —1. (We note that the chirality 
of these spirals for a given sign of the winding number is opposite of those in Ref. [17], 
because the reference condensate is on the outside rather than the inside of the ring.) 

In addition to generating the static trap, the DMD can also produce perturbations 
to the potential. However, because any individual mirror is binary (either on or off), 
potentials that require intermediate values require a form of grayscale control. We 
achieve such control by using Jarvis halftoning [28]. Because the point-spread function 
of our imaging system (in the plane of the atoms, ~ 6 pm 1/e^ full-width) is much 
larger than the DMD pixel size (in the plane of the atoms, ~ 0.5 pm), the potential at 
any given location is the convolution of the binary values of all the nearby pixels with 
the point spread function. For example. Fig. 2(b) shows a halftoned pattern written to 
the DMD for generating a sinusoidal perturbation of the form Vq sin(d) Convolution 
with the point-spread function generates the desired sinusoidal potential, as shown in 
Fig. 2(c). 

We empirically And that a perturbation of Vq ~ 0.4p applied for 2 ms is sufficient 
to excite the first phonon mode without perturbing the flow state of the ring. The 
speed with which our perturbation can be turned on and off is limited by the refresh 
rate of our DMD to 250 ps. To image the resulting density modulation in-situ, we 
use partial transfer absorption imaging (PTAI) [22], which is a type of minimally- 
destructive imaging [29]. In general, we transfer approximately 5% of the atoms into the 
imaging state. Despite the method being minimally destructive, none of the experiments 
contained in this work use repeated imaging on the same condensate. 

4. Results 

To detect rotation accurately in the ring, the behavior of the phonon mode in the 
absence of a persistent current must first be understood. To this end, we apply our 
perturbation to the inner ring and observe the resulting oscillation in the phonon 
mode, without a persistent current (£ = 0). Shown in Fig. 3 is the normalized ID 
density nxD(t,6) where nxD(t,6) is the ID density measured a time t after the 
perturbation was applied and no,iD(^) is the ID density separately measured with no 
perturbation applied. The density is determined via / n 2 r>{t, r, 6 ) rdr, where the 2D 
density n 2 D{t,r, 6 ) is determined from imaging followed by interpolated conversion into 
polar coordinates. The Bogoliubov mode oscillates with time in the ring, but is also 
damped (as discussed below). At 6* ~ ±90°, maxima and minima appear immediately 
after the perturbation. At these angles, clear oscillations are seen as a function of time. 

^ The halftoning is more evident if one zooms in on Fig. 2b. 
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Figure 3. Plot of the normalized ID density of the ring niD(t, 6>)/no,iD(^) (see 
text) vs. angle, for different times after the perturbation. The solid blue lines show 
the experimental data and the red dashed lines show the fit to the data. Each trace 
represents one shot of the experiment. For these data, there is no persistent current 
present (^ = 0). 


At the nodal points 0 0°, 180°, oscillations are expected to be absent in a perfectly 

uniform ring. 

We fit the data to the function sin(cjt + (fi) sm{9 + ^> 2 ) + c, where a, r, cj, 

(f>i^ 2 ^ and c are fit parameters. The average best fit value of the oscillation frequency 
is 79.3(3) Hz. We can estimate the frequency of this fundamental phonon mode as 
^0 = Cs/R^ where is the speed of sound. Using the speed of sound for a narrow 
channel ia/2m — 5.3(5) mm/s [30] yields (jJqI27t — 74(7) Hz, in good agreement 

with the average best fit value. The average best fit value of the decay constant is 
r = 25.6(1) ms, so that by 100 ms, the oscillation amplitude is diminished well below 
the noise. 

The decay of oscillation could be caused by a number of effects, including Landau 
damping and Beliaev damping, i.e., four-wave mixing. We expect Landau damping to 
thermalize our excitations by scattering from other, thermal, phonons [31]. In this case, 
we expect that the quality factor Q = ujt will be independent of the mode number 
q [24]. By contrast, damping via four-wave mixing, called Beliaev damping, will result 
in a strong dependence of Q on mode number [32, 24]. Landau damping should also 
depend on the temperature, while Belieav damping should depend on the density, but we 
have not explored these dependencies Nevertheless, we observe a Q for the g = ±1 

+ Furthermore, it is not clear if the phonon spectrum in one-dimension will fulfill momentum 
conservation in a four-wave mixing process. 
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Figure 4. Plot of the normalized ID density as a function of azimuthal angle and 
hold time after the perturbation. The three panels correspond to having different 
persistent currents in the ring: (a) ^ = — 1, (b) ^ = 0, and (c) ^ = 1. For each time and 
winding number, there is one experimental shot (i.e., there is no averaging). The green 
lines show the expected precession of the antinodes of the standing wave, according to 
Eq. 3. 

mode superposition of 12.8(6). We also investigated the decay of the q = ±2 mode, 
which has a Q of 13.8(4), agreeing within the uncertainties with the Q for the decay of 
the ±1. Therefore, Landau damping appears to be favored over Beliaev damping as the 
dominant damping mechanism. 

Having studied the relevant features of the phonon mode without a persistent 
current, we can now use Eq. 3 to predict the precession of the anti-nodes and nodes 
in the presence of a persistent current. The antinode initially at 0 = —90° reaches 
maximum density at approximately 10, 24, and 35 ms; the antinode initially at 0 = 90° 
reaches maximum density at approximately 5, 17, and 29 ms (see Fig. 4) *. At these 
times, the location of these antinodes is determined by Eq. 3. Shown in Fig. 4(a), (b), 
and (c) are the cases where the winding number in the ring prior to the perturbation 
was I = +1, 0, and —1, respectively. Given the radius of the ring, we expect that these 
peaks should move at rate of ih/mR^t = [^ x 1.25(7) deg/ms]L This rate is shown by 
the green lines in Fig. 4, and follows the peaks well. 

To detect a flow velocity with sufficient precision to distinguish between adjacent 
winding numbers using a single image, one needs to balance the observation time with 

* The maximum contrast of the standing wave does not occur at time t = 0 because of the details of 
how the perturbation is applied. 
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Angle (degrees) 


Figure 5. Histogram of the number of occurrences of measured precession angles of 
one of the antinodes of the standing wave at t = 35.1 ms. The red curve is a fit of 
three, independent Gaussians. These Gaussians correspond to three different winding 
numbers: t = —1, £ = 0, and £=1 (left to right). 


the signal to noise ratio. While the angular deflection grows linearly with time, the 
amplitude of the oscillation, and thus its signal to noise, decreases roughly exponentially 
with time. Empirically, we And that time t = 35 ms is the best compromise; here, the 
difference in deflection for A£ = 1 is expected to be ~ 45°. To test our ability to 
determine a given winding number with a single image, we lit the maximum of the 
density profile for 614 individual repetitions initialized with winding numbers £ = 0, — 1 
or 1. We bin the results to form a histogram, shown in Fig. 5. The histogram shows 
three clear peaks, which we fit to Gaussians. Each corresponds to one of the three 
winding numbers. Each has a slightly different width because of the compression or 
expansion of the phonon wave in the azimuthally varying density profile of the ring. 

We can use the Gaussian fits to predict the confidence with which we can assign a 
winding number. Given the strong overlap of adjacent Gaussians, we have a confidence 
of ~ 90% of identifying a ^ = — 1, and a confidence of ~ 95% of identifying an £ = 1 state. 
These confidences can be made better by attempting to make the ring more uniform to 
keep the form of the oscillation well defined. This might improve the fitting algorithm 
and decrease the noise. In addition, if the condensate were colder, the damping of the 
phonon mode should be reduced. This would allow for longer interrogation times and 
therefore more angular displacement between adjacent £ states. 

For this this method of measuring winding number to be minimally destructive, 
one must use a minimally destructive imaging method [29]. Previous experiments [33] 
show that such imaging methods do not perturb the flow state of BEGs. In addition, 
the perturbation one applies to the condensate must also be minimally disruptive. The 
observed fast decay helps to ensure this. Perhaps most important, the application of 
the perturbation must not change the winding number. To verify this, we measured the 
winding number that results from stirring at multiple different rotation rates using our 
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destructive interference technique. (This produced data similar to those in Ref. [3].) We 
then repeated this experiment, but after stirring, applied our sinusoidal perturbation, 
allowed 100 ms for the ring to settle, and measured the winding number again. The 
data show no significant change in the winding number as a result of the application of 
the perturbation. 

One might ask if going to higher modes, particularly q = ±2, would yield better 
sensitivity to rotation. The frequency of the oscillation increases as q, and because the 
quality factors of the modes are the same, the time at which one achieves a given signal- 
to-noise ratio scales as t ~ q~^. Because the precession angle scales as ~ t, these two 
factors cancel. However, one might be able to better determine the angular position of 
the maxima because the width of the oscillation peaks decreases as q~^. Further work 
is needed to verify these scalings. 

5. Conclusion 

We have demonstrated a minimally destructive in-situ method of measuring winding 
number in a ring-shaped Bose-Einstein condensate. This technique can be used as an 
excitation-based Sagnac interferometer [24] f. Because of its non-destructive behavior, 
this method can be applied to a variety of ring experiments, and may possibly be used 
to increase the sensitivity of atomtronic rotation sensors, by being able to repeat the 
measurement of winding number many times on a single condensate. 
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